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Abstract. Discrete strip-concave functions considered in this paper are, in fact, 
equivalent to an extension of Gelfand-Tsetlin patterns to the case when the pattern 
has a not necessarily triangular but convex configuration. They arise by releasing one 
of the three types of rhombus inequalities for discrete concave functions (or "hives") 
on a "convex part" of a triangular grid. The paper is devoted to a combinatorial 
study of certain polyhedra related to such functions or patterns, and results on faces, 
integer points and volumes of these polyhedra are presented. Also some relationships 
and applications are discussed. 

In particular, we characterize, in terms of valid inequalities, the polyhedral cone 
formed by the boundary values of discrete strip-concave functions on a grid having 
trapezoidal configuration. As a consequence of this result, necessary and sufficient 
conditions on a pair of vectors to be the shape and content of a semi-standard skew 
Young tableau are obtained. 

Keywords: Triangular grid, Gelfand-Tsetlin pattern, Discrete concave function, Younj 
tableau 



1 Introduction 

Let n G N. Consider a two-dimensional array X = (xjj)o<i<n,ai<j<6i of reals, where the 
index bounds a i? bi (depending on rows) are integers satisfying < 6j and: 

ao = 0, < <2i — do < <i2 — cti < . . . < a n — a n _i < 1, 

and 1 > 6i - b > b 2 - b x > . . . > b n - 6 n _i > 0. (1) 

We denote the set of pairs ij of indices in X by V and say that X has convex configuration. 
(This term is justified by the fact that V can be identified with the set of nodes of a convex 
triangular grid; see Remark 1 below. We visualize X so that (xoo, • • • ,^ob ) i s ^ ne topmost 
row and each triple Xy, x i+ ij, x i+ ij + i or Xy, a^j+i, x i+ xj+i is disposed so as to form an 
equilateral triangle. Then the array is shaped like a convex polygon, with 3 to 6 sides.) 
Two examples of such arrays are depicted in Fig. ^ 
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2 3 5 7 

2 4 5.5 4 , 1 6 9 11 

a: b: 

0354 -6 -1 356 

588 -8 -2 2567 



Figure 1: (a) a hexagonal array with n = 3, a = (0, 0, 0, 1), b = (2, 3, 3, 3); (b) a trapezoidal 
array with n = 3, a = (0, 0, 0, 0), b = (2, 3, 4, 5). 

Depending on the shape of the corresponding convex polygon, we may speak of hexag- 
onal configuration, pentagonal configuration, and etc. Although main results in this paper 
will be applicable to any of these, three special cases with a\ = . . . = a n = are of most 
interest for us: (a) hi = % for each % (giving a A-array); (b) b{ = i + m for each i (a 
array), see Fig. ^p; (c) b,i = m for each i (a [j- array), where m G N. In these cases we will 
also refer to an array as having triangular, trapezoidal, or parallelogram-wise configuration, 
respectively (usually ignoring other possible dispositions of triangle, trapezoid, or parallel- 
ogram). We say that X has size n in case (a), and (n,m) in cases (b),(c). Sometimes we 
will admit m = in case (b), regarding A-arrays as a degenerate case of ^-arrays. 

Let us associate with X the array dX = (&£#)o<i<ri,ai+i<j<&i of local differences dxij : = 
Xij — Xi,j-i, referring to dX as the row derivative of X. We deal with arrays X satisfying 
the following condition: for i = 1, . . . ,n and j = + 1, . . . , &j, 

dx^ > dxi-ij (when j < and dx^ij > ctejj+i (when j < bi). (2) 

The array dX obeying and having triangular configuration is said to be a Gelfand- 
Tsetlin pattern, and in this paper we apply the same name to dX with such a property 
when X has an arbitrary convex configuration as well. In this case we call X a strip-concave 
array, using an analogy with the corresponding functions explained in Remark 1 below. For 
example, both arrays in Fig. are strip-concave; their row derivatives are shown in Fig. [21 

2 1 5 2 

2 1.5 -1.5 , 5 3 2 
a: b: 

3 2 -1 5 4 2 1 

3 6 4 3 1 1 



Figure 2: Gelfand-Tsetlin pattern examples: (a) dX for X in Fig. EJi; (b) dX for X in 
Fig. Eh. 

One can identify the set of all arrays for V with the Euclidean space MX whose unit 
base vectors are indexed by the pairs ij G V. Let SCy denote the set of arrays X G M y 
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that satisfy property (J2J and the normalization condition x o = 0; imposing this condition 
leads to no loss of generality in what follows. Then SCy is a polyhedral cone in M y . 

Remark 1 . Let a, (3 be linearly independent vectors in R 2 . By a convex (triangular) grid 
we mean a finite planar graph G = (V, E) embedded in the plane so that each node of G 
is a point with integer coordinates in the basis (a, /3), each edge is the straight-line 
segment connecting a pair u, v of nodes with u — v G {a, /3, a + /?}, each bounded face is a 
triangle with three edges (a little triangle of G), and the union TZ of bounded faces covers 
all nodes and forms a convex polygon in the plane. A convex grid can be considered up to 
an afline transformation, and to agree with the above visualization of arrays, one should 
take the generating vectors as, e.g., a := (—1/2, —y/3/2) and (3 := (1, 0) and assume that 
(0,0) G V and > (0,0) for all G V. (The convex grids behind the arrays in 

Fig. ^ are exposed in Fig. El) A function x : V — > R determines an array X of convex 
configuration in a natural way: Xij := x(i,j). The arrays in SCy (considering V as the 
index set) are determined by the functions x having the following property: if / is the 
extension of x to TZ which is affinely linear on each bounded face of G, then / is a concave 
function within each region (strip) confined by the boundary of G and lines ia + WL/3 and 

(i — l)a + M./3, i — 1, 2, We call such a function x discrete strip-concave (by an analogy 

with discrete concave functions; see Remark 2 in the end of this section), and accordingly 
apply the adjective "strip-concave" to the arrays with property (J2J). 



a: 




Figure 3: (a) the grid for the array in Fig. ^i; (b) the grid for the array in Fig. ^p. 



Local differences on the "boundary" of X will be of most interest for us in this paper. 
These are represented by four tuples X x , X , v x (concerning the lower, upper, left and 
right boundaries, respectively) defined by 

\f := dx nj , j = 1, . . . , b n ; X f := dx 0j >, j' = 1, . . . , b ; 

(A vanishes when 6 — 0.) For example, the array X in Fig. has X x = (3, 0), A = (2, 1), 
H x = (2,-2,5) and v x = (1,0,4), and the array X in Fig. Cb has X x = (6,4,3,1,1), 
X X = (5, 2), fi x = (1, -7, -2) and v x = (4, -5, 1). 

Given A = (A a „+i, . . . , A 6n ), A = (Ai, . . . , X bo ) and \x,v G W 1 , define 

SC(X\X,n,u) := {X eSC v : (A x ,A X ,/i X ,z/ x ) = (X,X,fi,v)}. 
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This set, if nonempty, forms a bounded polyhedron (a polytope) in M. v in case of and 
£7-arrays. Indeed, (J2J) and £00 = imply 

Xij < /ii + . . . + Hi + Ai + . . . + Xj and x% > Hi + . . . + Hi + <?> (3) 

where q := A„_j + i + . . . + X n -i+j for ^Q-arrays, and q := Ai + . . . + Xj for ^/-arrays. (On 
the other hand, such a polyhedron V is unbounded when there is at least one interior entry 
and both left and right boundaries make a bend, i.e., < a n < n and < b n — b < n; 
in particular, if the hexagonal configuration takes place. One can check that adding any 
positive constant to all interior entries of an array IgP gives a point in V as well.) 

The first problem we deal with in this paper is to characterize the set By of all quadruples 
(A, A, Hi v ) (depending on V) such that SC(X \ A, /i, v) is nonempty. Two conditions on such 
quadruples are trivial. The first one comes up from the fact that (J2J) implies that A^ is 
weakly decreasing, i.e., X x +1 > . . . > X x , and similarly for A. The second one comes up 
by observing that 

\X X \ - |A X | + \/J, X \ ~\lS X \ = (X nbn - X nan ) - (Xob ~ X 00 ) + (x nan - X 00 ) - (x nK - X 0bo ) = 0, 

where for a tuple (vector) d = (d p , . . . , d q ), \d\ stands for ^2(di : % — p, . . . , q). 

To obtain the desired characterization, we need to introduce certain values depending 
on A, A. For k G Z + , define 

5 k (j) := max{0, \j_ k - Xj}, j = a n + 1, . . . , b n , and A fc := 5 k (a n + 1) + . . . + 6 k (b n ), 

letting by definition 5 k (j) := if j — k < or j — k > b . We refer to A k as the k-th deficit 
of A \ A. 

We shall explain later that the above problem is reduced to the case of trapezoidal 
configuration. Necessary and sufficient conditions on the corresponding quadruples for 
^Q-arrays are given in the following theorem. Hereinafter, for d = (d p , . . . , d q ) and I C 
{p, . . . , q}, d(I) denotes ^2(di : i G /), and for p < k < k' < q, d[k, k'} denotes d k + . . . + d k >. 

Theorem 1 For n G N and m G Z + , let X = (A 1; . . . , X n+m ) and X = (Ai, . . . , A m ) be 
weakly decreasing, and let /x, z/ G M n be such that \X\ — \X\ + \p,\ — \v\ = 0. Then a strip- 
concave £\-array X with (X x , X , fi X , v x ) = (A, A, //, v) exists if and only if the inequality 

A[l,|/|]+M/)-K/)-A|/|>0 (4) 

holds for each (including empty) subset I C {1, . . . , n}. Furthermore, if X, A, /i, v are integer 
and the polytope SC(X \ A,/i, u) is nonempty, then it contains an integer point. 

In particular, By forms a polyhedral cone (in M n+m x M m x M. n x IR n ) for V in question. 
Also (HJ) implies evident relations Aj > Xj (j = 1, . . . , m) and Xj < Xj- n (j = n + 1, . . . , n + 
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m), where the former is easily obtained by taking 7 = 0, and the latter by comparing 
|A| - |A| + \fi\ - \u\ = with © for / = {1, . . . , n}. 

Note that relation (jJJ involves a piece-wise linear term, namely, A|/|. One can replace 
each instance of (J3J) by a collection of 2 m linear inequalities, yielding an equivalent version 
of Theorem [T] It turns out that typically these inequalities determine facets of the cone By 
for £^-case; the precise list of facets of this cone is established in Section El (We shall see 
that the number of facets grows exponentially in n, m. On the other hand, to verify that a 
given quadruple (A, A, /i, v) belongs to By, it suffices to check validity of (JU) only for n + 1 
sets J: for k — 0, . . . , n, take I with |/| = k maximizing [y — //)(/).) 

For an arbitrary convex configuration, the problem with prescribed local differences 
A, A, /i, v is reduced to the trapezoidal follows. Since the polyhedron V := SC[X \ 

A, /x, z/) is described by a linear system formed by the inequalities in (j2J) and the corre- 
sponding equalities involving A, A, /i, z/, one can efficiently compute a number c G K + such 
that if V is nonempty, then there exists X G V with \xij\ < c/2 for all entries Xij. (For 
example, one can roughly take c equal to |V|'^' times the maximum absolute value a of 
the entries in A, A, fi, u, taking into account that the constraint matrix of the system has 
entries 0,1,-1. In fact, there is a bound c linear in a\V\; cf. (JSJ) for £7- arra y s -) Suppose 
a n and take the maximum p with a p = (then a« = i — p for p < i < n). Add to 

V the set A of pairs ij with 0<j < i — p < n — p, define A'- := c for j = 1, . . . , n — p, 
and define := /ij — c for i — p + 1, . . . , n. Symmetrically, if b n < b + n, we take the 
maximum q with b q = bo + q, add the set B of pairs ij with l<j — b n < i — q < n — q, 
define X'j := — c for j = 6 n + 1, . . . , b n + n — q, and define v\ := Vi — c for i = q + 1, . . . , n. 
Let A' coincide with A for the remaining entries, and similarly for v' . The resulting 

V := V U A U B gives a trapezoid (of size (n, &o))j and it is straightforward to verify that 
the set V := SC(X'/X, //, z/) (concerning V) is nonempty if and only if V is so, that the 
restriction of any X' eV toV belongs to V, and that X as above is extended in a natural 
way to an array in V' . 

Applying this reduction to the parallelogram- wise configuration of size (n,m), one can 
derive the following corollary from Theorem ^ 

Corollary 1 Let n,m G N, and let /i, v G M. n and weakly decreasing A, A G M. m satisfy 

y 

|A| — |A| + |/i| — \u\ = 0. Then a strip-concave o - array X with (X x , A , fi x , u x ) = (A, A,/x, z/) 
exists if and only if for each subset I C {1, . . . , n}, the inequality 

A[l, |/|] - X[m - |/| + 1, m] + /i(/) - z/(/) - A m > 

/ioWs /or |/| < m, and the inequality 

\X\ - | A| + - > 

/io/ds /or |/| > m. Furthermore, if A, A, /i, z/ are integer and the polytope SC(X \ A, /x, z/) zs 
nonempty, then it contains an integer point. 
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(To see this, observe that each entry v[ for the new right boundary tuple is equal to 
Vi — c, that //' = /i, and that A^- = — c for j = m + 1, . . . , m + n. The fact that A has all 
entries greater than — c implies that for k = 0, . . . ,n, each j with maxjm, k} < j < m + k 
contributes \j-k + c units to the new /c-deficit A' k (whereas S' k (j) = 5k(j) for j = 1, . . . , m 
and S' k (j) = for the remaining j's). Therefore, given I C {1, . . . ,n}, the new |/|-deficit 
becomes A|/| + A[m+1 — |/|, m] + |/|c whenever |/| < m, and |A|+mc whenever |/| > m. Also 
A'[l, |/|] = A[l, |/|] if |/| < m, and A'[l, |/|] = |A| - (|/| -m)c if \I\ > m. Now Corollary [His 
obtained from Theorem^ by substituting these relations, together with //(/) = /i(J) and 
v'{I) = — \I\c, into relation (j3J) (taken with primes).) 

A converse reduction, from [2r to EJ-c&se, is easily constructed as well, and Theorem^ 
follows from Corollary ^ In contrast, we cannot point out a "simple" reduction of Theo- 
rem^ to its special case with m = concerning A-arrays. (Nevertheless, a more intricate, 
though constructive, way of reducing does exist, as we explain in part D of Sectional In 
fact, this sort of reduction is behind our method of proof of Theorem ^ where the case 
m = is used as a base.) 

Another object of our study is the set of vertices of the polyhedron formed by strip- 
concave arrays X with convex configuration whose entries are fixed only on the lower, upper 
and left boundaries. More precisely, for A = (A an+i , . . . , Af, n ), A = (Ai, . . . , Aj, ) and /i 6 W 1 , 
define 

SC(X\X,fj) := {X e SC V : (\ x ,\ X ,fi X ) = (A,A,/i)}. 

(This polyhedron is bounded in case of A-, AV, or ^-configuration since the bounds on 
Xij indicated in (JHJ) remain valid in this case too.) We show the following. 

Theorem 2 For an arbitrary convex configuration and integer A, A, fi, the polyhedron SC(X\ 
A,//) is integral, i.e., each face of this polyhedron contains an integer point. 

Note that for arbitrary reals q±, . . . , q n , the transformation of an array X into the array 
X' with entries x'^ := Xy + qi preserves the row derivative. Such a transformation shifts a 
polyhedron SC(X\X, /i, u) into SC(\\\, /i', z/') with /4 := Hi + qi — qi-i and v[ := Vi + q% — qi~i 
(letting q := 0) and it maintains relation (0J). This implies that, without loss of generality, 
in Theorem^ one can consider only the quadruples of the form (A, A, n , i>) (where 0™ is 
the zero n-tuple). Similarly, one can restrict fi to be n in Theorem 121 as well. 

When dealing with A-configuration, for a triple (A, CP, z/), inequality (@J) turns into the 
majorization condition A[l, |/|] > Therefore, for a fixed A, the set {v : (X,0 n ,u) G 

By} forms a permutohedron, a polytope V formed by all vectors z G M. n with the same 
value \z\ such that for k = 1, ... ,n — 1, the sum of any k entries of z does not exceed 
a constant depending only on k. (The vertices of V are obtained by permuting entries 
of a fixed n-vector h; in our case, h = A.) It is known that for nonnegative integer A, z/, 
the majorization condition is necessary and sufficient for the existence of a semi-standard 
Young tableau with shape A and content z/, and that these tableaux one-to-one correspond to 
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the integer Gelfand-Tsetlin patterns respecting A, u; for a definition and a survey, see 
Theorem ^ (Corollary ^) shows that in case of ^-arrays (resp. £7 _arra y s ) an d A, A fixed, 
the analogous set {u : (A, A, CP, v) e By} forms a permutohedron in R n as well (but now 
the corresponding vertex generating vector h becomes less trivial to write down; it will be 
indicated in Section EJ). Each integer (generalized) Gelfand-Tsetlin pattern for nonnegative 
integer A, A, v determines a so-called semi-standard skew Young tableau with shape A\A and 
content v (cf. [H]), and our theorem (corollary) yields necessary and sufficient conditions 
for the existence of such tableaux. Figure 0] illustrates an instance of semi-standard skew 
Young tableau. 
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Figure 4: the semi-standard skew Young tableau corresponding to the pattern in Fig. Eb 
(here A = (6, 4, 3, 1, 1), A = (5, 2) and v = (3, 2, 3)). 

It should be noted that in case of A-configuration one can obtain the claim of Theorem|21 
by using a description for the generators of the Gelfand-Tsetlin patterns cone given in 
Berenstein and Kirillov pQ. 

Our method of proof of Theorem |21 is based on attracting a certain equivalent flow model 
and showing that the integer points in SC(\ \ A, n ) one-to-one correspond to the integer 
flows in a certain directed graph. In addition, we explain how to use the flow approach 
to easily show that Kostka coefficient K(X, u) (or K(X \ A, u)), as well as the intrinsic 
volume of SC(\, n , v) (resp. SC(\ \ A, n , z/)) in the nondegenerate case, preserves under 
a permutation of the entries of v. Here K(X, v) is the number of semi-standard Young 
tableaux with shape A and content v (which is equal to the number of integer points in 
SC(X, n , z/)), while K(\ \ A, u) concerns the corresponding skew tableaux. 

This paper is organized as follows. Theorems ^ and 121 are proved in Sections 121 and EJ re- 
spectively. Section El is devoted to a sharper version of Theorem [T] that precisely describes 
the set of linear inequalities determining facets of the cone By in ^Q-case (Theorem EJ). 
The concluding Section EJ discusses some additional aspects related to these theorems and 
demonstrates consequences from the proving method of Theorem EJ a combinatorial char- 
acterization of the vertices of polyhedra SC[\ \ A, //), the above-mentioned facts on integer 
points and volumes, and others. 

We conclude this section with two more remarks. 
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Remark 2 . Let us say that an array X (as in (jTfl) is (fully) concave if it satisfies (J2J) and 



Xij — Xi + ij > — for all 1 < i < n and a,i < j < b{. (5) 

This is equivalent to saying that the extension / of the function x on the nodes of the 
corresponding grid G (cf. Remark 1) is concave in the entire region 1Z. The functions 
x with such a property are often called discrete concave ones, and a series of interesting 
results on these have been obtained. Knutson, Tao and Woodward [Oj pointed out the 
precise list of facets of the cone BNDR n formed by all possible triples (A, /i, u) of n-tuples 
whose entries are the differences x(v ) — x(u) on boundary edges uv for a discrete concave 
function x on the triangular grid of size n, or a hive (equivalently: A, fi, v are the spectra 
of three Hermitian n x n matrices with zero sum). Also it is shown in [Hj that for each 
integer (A, /z, v) G BNDR n there exists an integer discrete concave function x as required 
for this triple. (A history of studying this cone and related topics are reviewed in see 
also (3|). Nontrivial constraints for BNDR n are expressed by Horn's inequalites. These are 
generalized to an arbitrary convex grid (see [S]), and relation (0J) in Theorem^is, in essense, 
equivalent to a special case of Horn's inequalites. We will briefly explain in Section El that 
Theorem Q can be derived from the above-mentioned results on discrete concave functions. 
At the same time, our direct proof of Theorem ^ is much simpler compared with the proofs 
of the corresponding theorems in jHJ E] . 

Remark 3 . The polyhedron integrality claimed in Theorem need not hold when the array 
entries are fixed on the whole boundary. More precisely, by a result due to De Loera and 
McAllister for any k G N, there exist A,/i, v G Z n and a triangular array X of size 
n, with n = O(k), such that A is a vertex of the polytope SC(X, fi, v) and some entry of 
X has denominator k. (Some ingredient from a construction in [3] is borrowed by [Zj to 
obtain an analogous result for fully concave triangular arrays in the case when the values 
are fixed only on two "sides".) Nevertheless, for /~\ -, £j- or A-configuration, at least one 
integer vertex in each nonempty polytope SC(\ \ A, fi, u) with A, A, /i, v integer does exist, 
as explained in the end of Section El 

2 Proof of Theorem H 

As explained in the Introduction, it suffices to consider the case /i = n . 

To show part "only if" in the theorem, we use induction on n. Case n = 1 is trivial, 
so assume n > 1. Let (A, A, 0™, v) e By (for V determined by n, m) and consider an array 
X G SC(X \ A, n , u) and a set I = {i(l), . . . , z(Jfe)} with 1 < i(l) < . . . < i(Jfe) < n. 

Define I' := I fl {1, . . . , n — 1} and X'j := 9x„_ij for j = 1, . . . , n + m — 1. Then 
Aj > X'j > X j+1 (by 0). By induction, 

A'[l, |/'|] - u(I') - Aj,,, > 0, (6) 
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where A' k , stands for the k'-th deficit for A', A, i.e., A' k , := S' k ,(l) + . . . + S k ,(n + m — 1), 
where S' k ,(j) := max{0, \j~k' — X'j}. Two cases are possible. 

Case 1. Let n I, i.e., I' = I. Since Sk(j) = max{0, \j~k — Aj}, S' k (j) = max{0, Xj-k — 
X'j} and Xj > A'-, we have 5 k (j) < 5' k (j), implying A fc < A' k . Now, using (JHJ), 

A[l, k] - v{I) -A k > A'[l, fc] - u{I) - A' k > 0, 

and © follows (with fi = n ). 

Case U. Let n £ I. Then |J'| = fc — 1. Summing up (JBJ) and the evident equality 
|A| — |A'| — v n — 0, we obtain 

A[l, k} + J2 ( A J " X U '■ j — k + 1, ■ ■ ■ ,n + m) - u(I) - A' k _ x > 0. (7) 

Note also that Xj + flfc(j) = max{Aj, Xj^k} (i n view of 5k(j) = max{0, Aj_fc — Xj}), and 
similarly \j_ t + ^^^j — 1) = max{A^_ 1 , Aj_jfc}. Since Xj < A^__ 1 , we have Xj + S k (j) < 
Xj-i + 8' k _ x {j - 1). Therefore, £(Aj - X'^ : j = k + 1, . . . , n + m) - A' fc _ 1 < -A k . This 
together with (JJJ) implies (jlj). 

Next we show part "if" in the theorem. We first consider case m = (i.e., A- 
configuration) ; in this case all deficits A& are zeros, which simplifies the consideration. We 
use induction on n; case n — 1 is trivial. Let n > 1 and let (jlj) hold for all /. In particular, 
Ai — v n > (by taking J := {n}). Also, subtracting inequality (jlj) with 7 = {1, . . . , n — 1} 
from the equality |A| — \v\ = 0, we obtain X n — v n < 0. Therefore, as A is weakly decreasing, 
there exists p G {1, . . . , n — 1} such that 

A p > v n and X p+1 < v n . (8) 
Assign the in — l)-tuple A' by the following rule: 

X'j := Xj for j = 1, . . . ,p - 1; := A i+i for j = p + 1, . . . , n - 1; 

and A^ := A p + X p+1 - v n . (9) 

Consider the triple (A', CP -1 , z/), where v' := {v\, . . . , z/„_i). We assert that 

A'[l, |/'|] > u(I') (10) 

holds for each /' C {1, . . . , n — 1}. Consider two cases, letting k := |/'|. 

(i) Let k < p. Then A'[l, k] = A[l, k), and (HDJ) follows from (jU) for / := I'. 

(ii) Let fc > p. Define I := I' U {n}. Then A'[l, jfc] = A[l, A; + 1] - v n (by ©), and we 
have (using (jlj)) 

A'[l, fe] - !/(/') = A[l, k+l)-u n - u(I') = A[l, |/|] - u{I) > 0. 
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Thus, HUH) holds for each I'. Also © and @ imply \ > X'j > X j+1 for j = 1, . . . , n - 1 
(in particular, A' is weakly decreasing), and Q together with |A| = \v\ implies |A'| = \v'\. 
By induction there exists a strip-concave A-array X' of size n— 1 with (X x> , fi x ' , u x ') = 
(A', CP -1 , z/). Assign := a;^- for < j < i < n — 1 and x n j := A[l, j] for 1 < j < n. The 
resulting array X of size n satisfies (j2J) and has the desired local differences on the "sides" , 
namely, (X x , \i x ,v x ) = (A,O n ,z/). Hence (A,O n ,z/) e £>y. Also when A, z/ are integer, the 
tuple A' defined by (JHJ) is integer as well, and the last claim in the theorem (for m — 0) 
follows by induction, as the integrality of X' implies that for X. 

It remains to prove part "if" when m > 0. Notice that the triple A, A, v can be considered 
up to adding a constant to all entries (which matches adding a constant to the array row 
derivative), so one may assume that A is nonnegative. Also, by compactness and scaling, 
w.l.o.g. one may assume that A, A, v are integer (this slightly simplifies technical details). 

We proceed by induction on m + |A|; case |A| = is trivial. Let (@J) hold for all /. In 
particular, Xj > Xj > Xj +n for j = 1, . . . , m. If X n+m = X m , we make a simple reduction to 
^-configuration of size (n, m — 1) by truncating the tuples A, A to A' := (Ai, . . . , A n+m _i) 
and A := (Ai, . . . , A m „i), respectively. (This maintains (jH), and if A is a required array 
of size (n, m — 1) for A', A ,u, then adding to X the elements Xi t i +m := Xi^+ m -\ + A m 
for i = 0, . . . ,n produces a required array of size (n,m) for A, A, v). A similar reduction 
(discarding Ai, Ai) is applied when Ai = Ai. 

Therefore, one may assume Ai > Ai and X m > X n+m . Then there are 1 < r < n + m — 1 
and 1 < s < m such that 

A r > Ai = . . . = X s > X s+ i and X s > X r+ i, (11) 

letting A m+ i := 0. Note that A r > A s+ i implies r < s + n and X s > A r+ i implies r > s. 
Define 

y . = f Aj - 1, j = r - s + 1, . . . , r, « 12 . 
j ' \ Xj, j = 1, . . . ,r - s, r + 1, . . . ,n + m; 



3 \ Xj, j = s + l,...,m. K J 

Then A', A are weakly decreasing and |A'| — |A | — \v\ =0. We assert that for any 
I Q {1) ■ ■ ■ , n} and k := |/|: 

X'[l,k]-u(I)-A' k >0, (14) 

denoting by A' k the fc-deficit for A', A , i.e., the sum of numbers S' k (j) := max{0, A-_ fc — A^} 
over j. To see this, first of all observe that S' k (J) = Sk(j) = if 1 < j < r (since A^ > Ai 
and X'j > J[, by dill),©,©)- Consider three cases. 

(a) Let k < r — s. Then for j — r + 1, . . . , n + m, we have A^- = Xj and X'j_ k = Aj_fc (in 
view ofj — k> s). Hence A' k = A k . Also A'[l, k) = A[l, k]. Then (JHJ) follows from (@J. 
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(b) Let r — s < k < r. Then S' k (j) — Sk(j) — 1 for j = r + 1, . . . , k + s (as 1 < j — k < s 
implies Xj_ k = Xj-k — 1 > A j — A^), and S k (j) = 6k(j) for j = k + s + 1, . . . , n + m. So 

= A fc - (k + s - r). Also A'[l, k] = A[l, k] — (k + s — r), and O follows. 

(c) Let r < k < n. Then ^(j) = 8k(j) — 1 for j = A; + 1, . . . , k + s, and S' k (j) = Sk(j) 
for j = k + s + 1, . . . , n + m. So A' k = A k - s. Also A'[l, k] = A[l, k] - s, and flHJ follows. 

Since |A'| < |A|, by induction the set SC(X' \ A ,0 n ,z/) is nonempty and contains an 
integer member X'. We transform X' into the desired array X for A, A, v as follows. Let 
a := \' r _ s+ i (=A r _ s+ i — 1). For i — 0, . . . ,n, define p(i) to be the maximum j such that 
dx'^ > a, letting by definition dx' i0 := oo. Then p(0) = 0, p(n) = r — s and p(i) < i for 

each i (as X' r _ s > X' r _ s+ i > \[ > dx' i i+1 ). For i — 0, . . . , n, define 



j = 0,...,p(i), 

a;y := ■{ ^+j-p(z), j = p(i) + 1, • • • ,p(i) + s, (15) 
x'ij + s, j = p(i) +s + l,...,i + m. 

Observe that A x = A, A = A and z/ x = (since x^j+m = ^ i+m + s for each i). Also 
A satisfies (J2J). To see the latter, let := — for all corresponding then 

G {0, 1}. Using the definition of a,p(0), . . . ,p(n), relation (fTo^) and the fact that X' is 
strip-concave, it is not difficult to conclude that €ij < ti-ij is possible only if j = p{i) = 
p{i — 1) + 1. In this case we have dx'^ > a + 1 > dx'^j, whence dxij > dxi-ij. Similarly, 
one can see that if ej_ij < e^+i, then j = p{i — 1) = p(i); in this case > dx it j + i 

follows from dx'^j > a + 1 > dx[^ +1 . This implies that X is strip-concave. 

This completes the proof of Theorem ^ 

Remark 4 . The A-array X recursively constructed in the second part of the proof is, in fact, 
a vertex of the polytope SC(X, n , v). This can be seen as follows. Given X' 6 SC(X, n , z/), 
let Q(X') be the set of all equalities of the form dx'^ = dx'^j or dx / i _ 1 j = dx' i:j+1 . A 
trivial observation is that X' is a vertex of SC(X, n , v) if and only if X' is determined by 
Q(X'), i.e., there is no other point X" in this polytope such that Q(X") I) Q(X'). In our 
case, the equalities as in Q (in the recursive process) give the corresponding equalities for 
dX; clearly the latter equalities determine X uniquely, so A is a vertex of SC(X,0 n ,i'). 
Moreover, if A, v are integer, then A is integer as well. This strengthens the last claim in 
the theorem for case m = 0. On the other hand, the construction of _arra y A in the third 
part of the proof does not guarantee that this A is a vertex of SC(X \ A, 0™, v). (Although 
an integer vertex in this polytope with A, A, v integer does exist, as explained in Sectional) 

Remark 5 . One can accelerate the process of constructing a required ,Cr arra y A in the third 
part of the proof. Given (not necessary integer) A, A, v, define p := Ai — max{A r+ i, A s+ i}, 
for r, s as in (JTTJ). When Ai > Ai and X m > X n+m , we can reduce the corresponding entries 
of A, A just by p (rather than by one), by setting A^- := Xj — p and A,- := Xj — p in the 
first lines of (fT2|) and (|TH|) . respectively (one shows that (jlj is maintained). Given an array 
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X' for A', A , v, we iteratively transform X' into an array for A, A, v. More precisely, at the 
first iteration, for a,p(0), . . . ,p(n) defined as above, we increase the entries x\a for ij as 
in the second and third lines of (fTHjl by e(j — p{i)) and by es, respectively, where e is the 
maximum value not exceeding p and such that the resulting array is still strip-concave (e 
is computed efficiently). If e < p, we apply a similar procedure (at the second iteration) 
to the updated X' and p := p — e, and so on. One shows that after 0(n 2 ) iterations we 
get p = 0, and that the final X' is the desired array X for A, A, v. Hence the number of 
operations in the whole process of finding a member of SC(X \ A, n , v) is polynomial in n. 
Such a transformation X' — > X is closely related to a rearrangement of flows (associated 
with strip-concave arrays) explained in part D of Section El 

3 Facets of the Cone By 

As mentioned in the Introduction, Theorem admits a reformulation in which the piece- 
wise linear constraints are replaced by linear ones. More precisely, one can see that for 
each / C {1, . . . , n}, inequality (HJ) is equivalent to the set of linear inequalities 

A[l, |/|] + A(J + |/|) - A( J) + //(/) - u{I) > 0, (16) 

where J ranges all subsets of {1, ... , m}, and for k G Z, J + k stands for the set {j + k : j G 
J}. In turns out that, as a rule, each of the latter inequalities is essential, i.e., determines 
a facet of the cone By. Note also that the description of this cone given in Theorem 
involves the "chamber inequalities" Xj > Xj + \ and Xj > Xj+%, so a priori such inequalities 
may determine facets as well. The precise list of facets of By is indicated in the following 
theorem. 

Theorem 3 For ^-configuration of size (n,m) with n > 1 and m > 0, inequality Mb\) 
determines a facet of By if and only if < \I\ + \ J\ <n + m and either (i) \I\ ^0,n (and 
J is arbitrary), or (ii) |/| = and \ J\ = 1, or (Hi) \I\ = n and \ J\ = m — 1. Furthermore, 
all these facets are different and By has no other facets if n = 1 or if n = 2 and m = 0. 
Otherwise the remaining facets are exactly those determined by the chamber inequalities 
Xj > Xj + \ for j = 1, . . . , n + m — 1, and Xj' > Xj> + i for f — 1, . . . , m — 1. In particular, 
By has (2 n — 2)2 m + n + 4m — 2 facets in case n > 2, and 2m facets in case n = 1. 

Proof. It is convenient to consider the reduced cone B* = B* nm formed by the vectors 
(A, A, v) G W n+rn x R m x R n such that {X,X,0 n ,u) belongs to the cone B := By. Since 
(A, A, /i, u) G B if and only if (A, A, n , v — p) G B (as explained in the Introduction), the 
cone B* lies in the hyperplane TC defined by |A| — |A| — \v\ = and is described by the 
inequalities 

A[l, |/|] + A( J + |/|) - A(J) - > (17) 
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for all J C {1, . . . , n} and JC {1, . . . , m}, and by the chamber inequalities. Moreover, there 
is a natural bijection between the facets of B and B* , namely: for any /, J, inequality (fTo^ 
determines a facet of B if and only if (fTTj) determines a facet of £>*, and similarly for the 
chamber inequalities. 

In light of these observations, our goal is to characterize those of the instances of (fTTjl 
and of the chamber inequalities that determine the facets of B*. We will argue in terms of 
the grid G = (V,E) (defined in Remark 1 in the Introduction). Let E* = E* m be the set 
of "horizontal" edges := — 1), of G. An edge with i = n (i = 0) is also 

denoted by ej (resp. e^); we refer to the sets L := {ei, . . . , e n+m } and L := {ei, . . . ,e n+m } as 
the /operand upper boundaries, respectively. Besides, we will deal with the edges r\, . . . , r n 
on the right boundary R, where Ti connects the vertices (i — 1, i + m — 1) and (i, i + m). 
We associate with an edge e G E* the function \ e ° n -S* taking value 1 on e, and on the 
remaining edges. 

Let h R. The border of ft, is defined to be the function (T = a' 1 onlULUi? 

coinciding with h on L U L and taking the value {h{e^i) + . . . + h(e iji+m )) — (h(ei-i t x) + 
. . . + hfe-ij+m-x)) on ri, i — 1, ... ,n; we also use vector notation, identifying a with the 

corresponding triple (\ h , X , u h ), where := cr(ej), := cr(ej) and := cr(rj). Also h is 
identified with the corresponding array dX of row derivatives (i.e., dxij = /i(e^) for all ij), 
and we say that h is a Gelfand-Tsetlin pattern, or, briefly, a GT-pattern, if <9X is such. In 
other words, h is a GT-pattern if /i(e) > h(e') for all pairs (e, e') of the form (e^, e^ij) or 
(e^, ej + i J+ i); we denote the set of these pairs by II. 

In what follows, when proving that one or another instance C of valid inequalities is 
facet-determining, we will try to construct 2n + 2m — 2 GT-patterns for G such that their 
borders are linear independent and each of them turns L into equality. First of all we show 
that the cone B* has co-dimension one (recall that B* is contained in the hyperplane 7i). 

Claim. B* nm contains 2n + 2m — 1 linearly independent vectors (assuming n > 1 and 
m > 0). 

Proof. Consider the 2n + 2m — 1 functions hi, ... , h n+m , hi, ... , h m , qi, . . . , q n ~i on E* 
defined by 

hj : = X ej ; hj'= X* 3 ; ft := X^ t+m ■ 
Observe that the border of each hj (hj) takes value 1 on the edge ej (resp. ej), and on 
the other edges in L U L U R, while the border of takes value 1 on rj_i, —1 on rj, and 
otherwise. Therefore, these 2n + 2m — 1 borders are linearly independent. However, the 
above functions, except for hi, are not GT-patterns (e.g., ^2(^-1,1) = < 1 = ^2(62))- By 
this reason, add ch to each of these functions, where c is a large positive number and h is 
the function on E* defined by 

h(eij) := z - 2j. 

One can see that h is a GT-pattern; moreover, h(e) — h(e') = 1 for each pair (e, e') G II. 
This implies that the updated functions are already GT-patterns. Also their borders remain 
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linearly independent (as c is large). So we have 2n + 2m — 1 linearly independent vectors 
in £>*, as required. ■ 

Next consider / C {1, . . . ,n} and JC {1, . . . , to} with < |/| + | J| < n + m. (When 
|/| = n and | J\ = m, (|T7jl becomes |A| — |A| — \u\ > 0, which is not essential.) We examine 
the three cases of I indicated in the theorem. 

Case 1. Let < |/| < n. We have to show that (j!7|) determines a facet for any J. 
Our construction of 2n + 2m — 2 linearly independent vectors that belong to B* and attain 
equality in (fT7j) is based on a certain partition of E* into subsets Ai, . . . , A n+m (depending 
on /). Each Aj is defined to be a minimal set satisfying the following conditions: 

(18) (i) Aj contains tj\ 

(ii) if Aj contains with ^ i G" / and p < i + m, then ej„i iP G Aj] 

(iii) if A,- contains e ip with j G / and p > 1, then ej_i iP _i G Aj. 

These sets are defined uniquely and do give a partition P of E*. (For example: if n — 3, 
m = 2 and / = {2}, then A 1 = {e 1 ,e 2 i}, A 2 = {e 2 ,e 2 2,e n ,e 1 }, A 3 = {e 3 , e 23 , e 12 , e 2 }, 
At = {e 4 , e 2 4, ei 3 }, A 5 = {e 5 }.) Also for j — 1, . . . , m, both edges e,, ej+\j\ are contained in 
the same set in P, namely, in Aj + \i\. 

Define the function h on E* by 



The following properties of h will be important for us: 

(20) (i) h is a GT-pattern, and h(e) > h(e') holds for each pair (e, e') G II such that e, e' 
belong to different sets in P; 

(ii) (fTTjl holds with equality for the border (A, A, u) of h. 

Indeed, (i) follows from (fTTI|) because if Ci V G Aj then each of ej_i iP and ej + i jP+ i (if any) 
belongs to Aj or A, +1 , by (fTHj) in the construction of P. To see (ii), let £j := n + m — j for 
j = l,...,n+TO. Then A[l, |/|] = 4 + and A(J+|I|) = ^(^ : j G J+\I\) = X(J). 

Observe that for each i G / and p = 2, . . . , i + m, both edges ejp, ej_i iP _i belong to the same 
set in P, and the edge e^x belongs to the set A^(j), where <i(i) is the number of elements of 
I greater than or equal to i. This implies 



whence u(I) = £(£ d (i) : i & I) = h + ■ ■ ■ +% Therefore, A[l, |/|] = u(I) and X{J+\I\) = 
A(J), yielding® (ii). 



/i(e) := n + m — j 



for e G Aj, j — 1, . . . , n + to. 



(19) 
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Next we construct 2n + 2m — 2 functions on E*, not necessarily GT-patterns, such 
that their borders are linearly independent and attain equality in (|T7j) . To this aim, we 
form two auxilliary £A,-grids G' = (V',E f ) and G" = (V",E"), the former having size 
{n' := \I\,m' := \J\) and the latter having size (n" := n — n',m" : = m — m'). For 
convenience, we use notation with primes (double primes) for edges and their sets in G' 
(resp. G"). Since < |/| < n, we have n', n" > (while m', m" > 0). By the Claim applied 
to G', there exist GT-patterns ai, . . . , a2 n '+2m>-i '■ E'* — > K (concerning G") whose borders 
are linearly independent. Similarly, there exist GT-patterns b\, . . . , &2n"+2m"-i : E * — > M 
(concerning G") whose borders are linearly independent. 

These patterns are transformed ("lifted") into functions on E* by use of special maps 
U\ : Ei — ► i?'* and uj 2 : £2 — ► -E 1 "*, where i?i (£'2) is the union of the sets Aj for j G Ji := 
{1, . . . , I J|} U ( J + |/|) (resp. for j G J 2 := {l,...,n + m} \ Ji). The map uii is defined so 
as to satisfy the following condition: 

(21) for ey, e P9 G Ei, e^-, := Wi(ey) and e' plq/ = uJi{e pq ), 

(i) if i = p and j < q, then z' = p' and j' < q'\ 

(ii) whenever e^, e P9 belong to the same set A d (d G J\): if p = i — 1 and g = j (and 
therefore, z G" /, by p8jl ) then i'j' = p'q'\ and if p — i — 1 and q — j — 1 (and 
therefore, i E I) then p' = i' — 1 and q' — j' — 1. 

In its turn, is defined so as to satisfy: 

(22) for aj, e pq G E 2 , e" f>jll := ^(ey) and e£,y, = ^2(e pg ), 

(i) if i = p and j < q, then i" = p" and j" < g"; 

(ii) whenever ey, e pg belong to the same set A rf (cf G J2): if p = i — 1 and q = j, then 
p" = i" — 1 and g" = j"; and if p = i — 1 and q = j — 1, then = 

One can check that both Ui, u)i are well-defined and unique. Also uoi establishes one-to- 
one correspondence between the sets Li := {ej : j G Ji} and L' := {e[, . . . , e^, +m ,}, as well 
as between the sets Li := {ej : j G J} and L := {e^, . . . ,e' m ,}. Similarly, uj 2 establishes 
one-to-one correspondence between L 2 := {tj : j G J 2 } and L" := {e", . . . , e^// +m «}, and 
between L 2 := L\ Li and L ' := {e'/, . . . , e^«}. 

Using and u> 2 , the above-mentioned functions a s and b t are transformed into functions 
on E* in a natural way, as follows. For s — 1, . . . , n' + m! — 1, define 




for e G 



for e G 



'2- 



(23) 



and for £ = 1, . . . , 2n" + m" — 1, define 




for e G 



'2- 



(24) 
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Let a s (( t ) be the border of g s (resp. h t ). The values of a s on the edges in L\ U L\ 
coincide with the values of the border a s of a s on the corresponding edges in V U L , and 
a s takes zero values on the remaining edges in L U L. Also (|21|) and (|23|) show that for 
i = 0, . . . , n, the sums 5^,1) + . . . + &,(ei ii+m ) and a s (e-, -J + . . . + a^e-,^^^) are equal, 
where i' is the number of elements of I smaller than or equal to i. This implies that the 
value of a s on the edge of the right boundary of G is zero if i g" I, and equals the 
value of a s on the edge r[, of the right boundary of G' if % G /. Similar correspondences 
take place for ( t and the border of b t (regarding the sets L 2 , L 2 , L", L and replacing i' by 
i»:=\{0,...,i}\I\). 

So we can conclude that these 2n + 2m — 2 vectors a s and are linear independent. 
Also the equality |A'| — |A | — \v\ — for each a s = (A', A , z/') implies A(Ji) — A(J) —v{I) = 
for cx s = (A, A, z/), while the latter equality holds automatically for each ( t = (A, A, u) since 
is zero within L\, L\ and {ji : i £ J}. Thus, each of the obtained borders turns (j!7|) into 
equality. 

Since the functions g s and /i 4 need not be GT-patterns, we add ch to each of them, 
where h is defined in (II 9|) and c is a large positive number. Then the new functions are 
GT-patterns and their borders are as required. (The former relies on property (|2T)|) and 
the relationships between g s ,(t and a s ,b t : observe that for (e, e) G II, if both e, e are in 
the same set of the partition P, then ^ s (e) > ^(e), h t (e) > h t (e) and /i(e) = h(e), while if 
they are in different sets, then h(e) > h(e).) This completes the proof for case (i) in the 
theorem. 

Case 2. Let 7 = 0. Then (fTTj) becomes A(J) — A(J) > 0, which is the sum of valid 
inequalities Xj — Xj > over j G J. So we have to examine only the latter inequalities 
for j = 1, . . . ,m (each concerning the case \J\ — 1). To show that each A^ — Xj > is 
facet-determining, we argue as in Case 1 and apply the reduction to the grids G f , G" for 
J = {j}. The grid G" has size (n, m — 1) and generates 2n + 2(m — 1) — 1 = In + 2m — 3 
linearly independent vectors in the intersection of B* and the hyperplane Xj — Xj = 0. 
Moreover, these vectors (A, A, z/) satisfy Xj = Xj = 0. One more vector comes up from the 
degenerate grid G' (having size (0,1)); it has entries Xj = Xj = 1, and otherwise. (This 
is the border of the GT-pattern taking value 1 on the edges e 0j -, eij, . . . , e n j, and on the 
remaining edges in E* .) This yields (ii) in the theorem. 

Case 3. Let \I\ = n. This case is symmetric to Case 2. Inequality ()17|) becomes 
A[l,n] + A(J + n) — A(J) — \v\ > 0. This is equivalent to — A(J + n) + A(J) > 0, where 
J := {1, . . . ,m}\ J (assuming |A| — |A| — \v\ = 0). The latter is the sum of valid inequalities 
—Xj +n + A j > over j E J (each being equivalent to (fTTj) with an (m — l)-element set as 
J). To show that for each j = 1, . . . , m, the inequality — Xj +n + Xj > determines a facet of 

we apply the reduction to the grids G',G" for J = {1, . . . , m} \ {j}. Then 2n + 2m — 3 
vectors come up from G', and one vector from G" . (Another method: use the symmetry 
of B* defined by (A, A, z/) — > — (A„ +m , . . . , Ai, X m , . . . , X%, v\, . . . , z/ n ), which reduces case (iii) 
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to case (ii) in the theorem.) 

It remains to examine the chamber inequalities. When n = 1, each inequality Xj > A J+ i 
is not essential, as it is the sum of Xj > Xj and Xj > A J+ i (which are the instances of (fTTj) 
with J := {j} and with J := {1, . . . , m} — {j}), and similarly for the inequalities Xj' > Xf+i- 
When n = 2 and m = 0, the unique chamber inequality Ai > A2 is not essential as well, as 
it follows from the inequality Ai — V\ > (i.e., (fT7|) for / := {1}), the inequality Ai — v 2 > 
(i.e., (fT7j) for / := {2}), and the equality Ai + A2 — V\ — V2 = 0. 

Now consider the case n > 2 and n + m > 3. We show that for j = 1, . . . , n + m — 1, 
Xj > Aj + i determines a facet of B*, as follows. Take the following vectors: 
the border of x ej for j = 1, . . . , j — 1, j + 2, . . . , n + m; 
the border of x £j for j = 1, . . . , m; 

the border of x eip f° r i = 1, ■ ■ ■ , w — 1 and one p such that 1 < p < i + m and 
(i,p) =£(n- 1, j); 

the border £ of the all unit function on E*. 

(An edge e n _i jP different from e n _ij exists since n + m > 3.) These 2n + 2m — 2 vectors 
are linear independent and satisfy the equality Xj = Xj + ±. We modify each of them by 
adding ca, where c is a large positive number and a is the border of a GT-pattern h on 
E* such that h(ej) = h(ej + i) = h(e n -ij) and h(e) > h(e') for all pairs (e, e') G II except 
for (ej,e„_ij) and (e„_ij, e^+i). (To construct such a pattern is easy.) The new vectors 
become the borders of GT-patterns and are as required. The fact that each inequality 
Xj > Xj + \ (when n, m > 2) determines a facet is proved in a similar way. 

We leave it for the reader to verify that all facets of £>* appeared throughout the above 
proof are indeed different. (In other words, for the facet-determining inequalities claimed 
in the theorem, their incidence (0, ±l)-vectors have the property that no distinct vectors 
£, £' satisfy £ = a(,' + j39, where a G R+, (3 G M, and 9 is the incidence vector of the equality 
|A| - |A| - \v\ = 0.) 

This completes the proof of Theorem El 

4 Proof of Theorem 2 

First of all we observe that the generic case of convex configuration in this theorem is 
reduced to the case of ^-configuration. Indeed, given A, A,/z for V as in ((H), there exists 
a (sufficiently large) positive integer c such that each face of SC(X \ A,/i) contains a face 
of the polyhedron V formed by the arrays X G SC(X \ X, fx) with \dxij\ < c for all entries 
dxij of dX. Let m := b and extend A to {n + m)-tuple A' by setting X[ :=...:= X an := c, 
K(n)+i ■= ■ ■ ■ ■= X' n+m := -c and A^- := A^ for j = a n + 1, . . . , b n . Accordingly, set li[ := ^ 
if ai = 0, and \i\ := fii — c if a« > 0. Then the restriction map X' X'\ y gives a bijection 

between the ^.-arrays X' with X x = A', A = A, fi x = fi' and the arrays in V (cf. 
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explanations in the Introduction). This implies that SC(X \ A, (i) is integral if SC(X' \ A, //) 
is such. 

In the rest of the proof we deal with ^-configuration of size (n, m). As before, we may 
assume fi = 0™. Also one may assume that A is nonnegative (cf. reasonings in the previous 
section). For brevity we denote the polytope SC(X \ A, n ) by SC(X \ A). Theorem 2 will 
be proved by constructing a bijection between the vertices of SC(X \ A) and certain forests 
in the grid G (defined in Remark 1 in the Introduction). Establishing this correspondence, 
we admit A and A to be real-valued. 

The node set V of G is naturally partitioned into subsets (horizontal layers) Li = {(i, 0), 
. . . , (i, i + m)}, i — 0, . . . , n. Extract the edges connecting neighbouring layers and orient 
them from the top to the bottom. Formally: let A be the set of pairs e°- := ((i,j), + 
and e]j := ((z, j), {i + 1, j + 1)) of nodes of G, for i — 0, . . . , n — 1, j — 0, . . . , i + to. Then 
H n ,m '■= H := (V, A) is an acyclic digraph in which any maximal (directed) path begins at 
a node of the "topmost" layer L and ends at a node of the "bottommost" layer L n . 

We say that a function g : A :— > R + is a (A, X)- admissible flow in H if 

{0, % = l,...,n- 1, j = 0, ...,z + m, 

Xj-Xj+i, i = n, j = 0,...,n + m, (25) 

A i+ i - Xj, i = 0, j = 0, . . . ,m. 

Here div g (v) (v G V) stands for the value J2 e =( u ,v)eA d( e ) -J2 e =(v,u)eA #( e )' and we formally 
extend A and A by setting Ao := Ao := Ai and X n+m+ i := A m+ i := 0. In particular, 
g(^n-i o) = anc ^ ^(^-i n+m-i) = ^n+m- The set JF(A \ A) of (A, A)-admissible flows forms 
a polytope in M' A L 

Claim. For any X e 5C(A \ A) tiere exists a (A, X) -admissible How g = j(X) satisfying 
1(f)- t i,j ~ 9Xi+ di +1 ' i-0,... ; n-l ; j-0 : ...A-m : (26) 

letting dx i0 := Ai and dxa +m+ i := 0. Moreover, 7 is a bijective mapping of SC(X \ A) to 
^(A\A). 

(Figure El illustrates the flow determined by the array X with dX as in Fig. |2|d; here 
the flow is integer and its value on an edge is indicated by the number of lines connecting 
the ends of this edge.) 

Proof. Let X G SC(X \ A) and let g be defined by f!26|) . Then for each node v = (n,j) 
with j — 0, . . . , n + m, 

div s (t)) = giel^^) + g(el_ hj ) = (dx nj - dx n - ltj ) + (dx n -ij - dx nd+1 ) = Xj - X j+1 , 

letting g(e) := if the edge e is void (e.g., for e = e\_ x Similarly, div g (t>) = A J+ i — Aj 
for each node v = (0, j), j — 0, . . . ,m. And for each node v = (i,j) with 1 < i < n — 1 and 
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Figure 5: the flow corresponding to dX in Fig. |2fc>. 

< j < i + m, one has 
div g {v) = g{e\_ ld _ x ) + g(^-u) - g{e%) - g{e\ 3 ) 

= (dxij - dxi-ij) + (dxi-ij - dxij+i) - (dxy - dx i+hj+1 ) - (dx i+1J+1 - dx iJ+1 ) = 0. 

Also the function g is nonnegative, as is seen by comparing ()26|) and (J2J). Thus, g is a 
(A, A)-admissible flow. 

Conversely, let g be a (A, A) -admissible flow in H. Assign numbers dx^ recursively by 
the following rule: 

dx nj := Xj, j = 1, . . . ,n + m, 

dx^ := dx i+1J - giejj^), i = 1, . . . , n - 1, j = 1, . . . , i + m. 

This gives the ^,-array X of size (n, m). Reversing the argument above, one can check 
validity of ([26)1 . This and the nonnegativity of g imply that X is strip-concave and satisfies 
X x = X and A = A. Then X £ SC(X \ A), and the claim follows. ■ 

Thus, 7 is a linear operator (in view of ([26))) and 7 gives a one-to-one correspondence 
between the points in the polytopes SC(X \ A) and JF(A \ A). Therefore, 7 establishes a 
one-to-one correspondence between the vertices of these polytopes. 

Next we characterize the vertices of JF(A \ A). To this aim, we distinguish, in the 
bottommost layer L n , the set L(X) of nodes (n,j) (1 < j < n + m) such that Xj > Xj + i, 
and in the topmost layer L , the subset L(X) of nodes (0, j) (0 < j < m) such that 
Xj > Aj + i. Given a flow (7 G JF(A \ A), let H(g) denote the subgraph of H induced by the 
set of edges e with g(e) > 0. From (j2^jl it follows that H(g) contains L(X) and L(X) and that 
each node of H(g) lies on a path from L(X) to L(A). Suppose there are two different paths 
P, P' in H(g) having the same beginning and the same end. Choose e > not exceeding 
the minimal value of g on the paths P and P' . Then the functions g' := g + ex P — €X P ' 
and g" := g — t\ F + e% P ' are nonnegative and satisfy pHjl . where x Q G {0, 1} A is the 
characteristic function of the edge set of a path Q. So g is expressed as the half-sum of two 
different (A, A)-admissible flows g',g", and therefore, g cannot be a vertex of JF(A \ A). 

On the other hand, let for any two nodes y and z, H(g) contain at most one path from 
y to z, i.e., H(g) is a (directed) forest with the set L(X) of zero indegree nodes (roots) 
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and the set L(X) of zero outdegree ones (leaves). Then g is the only (A, A)-admissible flow 
taking zero values on all edges outside H(g), i.e., g is determined by H(g). Indeed, one can 
see that for each edge e = (u, v) of H(g), g(e) is equal to 

- A i+ i : (n,j) e V(Q)) - £(X, - X j+1 : (0,j) e V(Q)), (27) 

where Q is the connected component of H (g) \ {e} that contains the node v, denoting by 
V(Q) the node set of Q. This implies that g is a vertex of JF(A \ A). Moreover, g is integer 
if A, A are integer, and Theorem |21 follows. ■ 

Arguing as in the above proof, one can associate the vertices of SC(X \ A) with certain 
subgraphs of H, as follows. 

Corollary 2 In case of /~\- configuration of size (n,m), each vertex of SC(X\X) one-to-one 
corresponds to a forest H' in H n ^ m having L(X) as the set of roots and L(X) as the set of 
leaves and satisfying the following condition: for each component Q of H' , the value in (fF7| ) 
is zero, and for each edge e = (u, v) and the component Q of H'\{e} containing v, the value 
in \2T\) is positive. Therefore, in case m = 0, the vertices of SC(X) one-to-one correspond 
to the rooted trees in H n (:— H n ) with root (0,0) and set of leaves L(X). 

Remark 6 . The flows introduced in the proof of Theorem 2 give an alternative way to rep- 
resent the Gelfand-Tsetlin patterns (or the strip-concave arrays), and Corollary |2] suggests 
a way to compute or estimate the number of vertices of the polytope SC(X \ A,/i) in case 
of /2r con fig ura tion (or A-configuration). One can check that the reasonings in the proof 
of Theorem |21 and the corresponding corollary are applicable to o- con fig ura ti° n as we U 
(with H n ^ m arising from the corresponding parallelogram- wise grid of size (n, m)). 

5 Concluding Remarks 

In this section we outline (in parts A-D) more applications of the flow approach developed 
in the proof of Theorem 121 Here, unless explicitly said otherwise, we consider the case of 
^-configuration of size (n, m). (Note that the exposed properties remain valid if we deal 
with ^/-configuration.) 

A. Let V = Vn t m be the set of paths in the graph H = H n)Tn beginning at a node of the 
layer L and ending at a node of L n \{(n, 0)}. Associate with a path P 6? the 0" arra y Y p 
with the entries y a = . . . = = 1 and y^+i = ... = yi, i+m = for i = 0, . . . , n, where 
(i,p(i)) is a node of P. Considering the case of triangular arrays, Berenstein and Kirillov [l 3 
noticed that the set of arrays Y p (P e V n fi) constitutes a minimal list of generators of 
the cone of nonnegative Gelfand-Tsetlin patterns of size n — 1. A similar property takes 
place for ^-patterns (or ^-patterns) and can be easily shown by use of flows. More 
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precisely, for a strip-concave £2r arra y X with dX > 0, take the flow g = 7(X) defined 
by (|26|). Then g is represented as a nonnegative linear combination a\X Pl + • • • + &nX Pn , 
where Pi, ... , Pn G P. One can check that dX = a\Y Pl + . . . + o>nY Pn , as required (the 
minimality of {Y p : P G P} is obvious). 

B. One can establish some invariants for polytopes SC(X \ A, 0", v) when the entries of 
v are permuted. Consider an array X G SC(X \ A, n , v) and the flow g = j(X) as in (j2Sj) - 
For z = 1, . . . ,n, we have ^tT^ij - ZljtT" 1 ^i-i.i = x i,i+m ~ Xj_i ji+m _i = z^. Also 
c^Xjj — dxi-ij = g{e\-ij-\) for j = 1, . . . , z + m (see Section |U for the definition of edges 
e°j-/ and e^-,; as before, (9xj_i 5 j +m := 0). Comparing these relations, we conclude that 

n = ^(e-_i,o) + • • • + ^(ej_i, i+m _i) for i = 1, . . . , n. (28) 

Choose i G {1, . . . ,n — 1} and consider the subgraph H % of H induced by the edges 
connecting the layers Lj_i, Lj or the layers Lj, For j = 0, . . . , i + m — 1, the nodes 

{i — 1, j) and {i + 1, j + 1) are connected by two paths, namely, by path Zj with the edges 
e i-iji e lj an d by path Zj with the edges ej-ij, e i,j+i- Let us call such a path Z with edges 
e, e' a zigzag and define its capacity to be g(Z) := min{g(e), ^(e')}. The zigzag swapping 
operation modifies g within H l by swapping the capacities simultaneously for each pair 
Zj, Z'y More precisely, for j = 0, . . . , i + m — 1, assign 

//a ._ f S'(e) - flr(Z/) + g{Z' j ) for each edge e of Zj, 

\ d{ e ) ~ 9(Z'j) + g(Zj) for each edge e of Zj, 

and g'(e) := g(e) for the remaining edges of H. Obviously, g' is again a (A, A)-admissible 
flow. (For example, such an operation applied to the flow in Fig. El results in the flow 
illustrated in Fig. EJ) 




Figure 6: the flow obtained by the zigzag swapping operation applied to the flow in Fig. El 
at layer 2. 

Let o~i(X) denote the array 7 -1 (g') and let v' be the n-tuple of local differences on the 
"right boundary" of o"j(X). Using (|2~%|). one can check that the zigzag swapping operation 
swaps Vi and z/ i+1 , i.e., v[ = u i+ i, u' i+l = z/j and v' — v v for p i, i + 1. Moreover, applying 
the zigzag swapping operation (with the same i) to g' returns g. 

Thus, for each i, o~i is a continuous bijective mapping of <SC(A\A, n , v) to SC(X\X, n , v') 
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(and of is the identity on SC(X \ A)) 1 . Moreover, for k G N, if g is ^-integer, so is g'. 
Therefore, Oj gives a bijection on the ^-integer points in these polytopes for any k. As a 
consequence (for k = 1), the following known property is obtained: if A, A, v are integer and 
if v' is an arbitrary permutation of u, then Kostka coefficients K(X \ A, v) and K(X \ A, i/) 
are equal. 

C. Let A, A, v be rational- valued and let v' be a permutation of v. Let Vo denote the set 
of boundary index pairs in V (or the boundary nodes in the grid G) . The fact that each 
map a", is continuous and bijective implies that the polytopes SC := SC(X \ X,O n ,u) and 
SC' := <SC(A\A, n , v') have the same dimension (which typically equals |V\Vo|). Consider 
the \V \ Vo (-dimensional affine subspaces S and S' containing the polytopes SC and SC' , 
respectively, which are obtained by imposing the corresponding equalities on the values on 
Vo- Since S and S' are parallel, there is k! G N such that for any multiple k of k', the lattice 
of -^-integer points in 5" is obtained by a parallel translation of a similar lattice in 5*. So 
the density of ^-integer points in S and S' (measured by the number of such points in a 
unit ball with center at a point of the lattice) is the same. Also the numbers of ^-integer 
points in the polytopes in question are equal. Thus, when k tends to infinity, we obtain 
equality for the corresponding volumes and can conclude with the following. 

Proposition 1 Given (real-valued) A, A, v , let v 1 be a 'permutation of v. Then the polytopes 
SC(X \ A, n , v) and SC(X \ A, n , z/) have the same \V \ V \- dimensional volume. 

It should be noted that, although o"j (being a piece-wise linear operator) brings integer 
points into integer ones, it need not do so for polytope vertices, even for polytopes SC(X\X). 
Indeed, in case m = 0, take a rooted tree T in H n (with root (0,0) and the leaves in L n ) 
such that for some i,j, the subgraph TdH" 1 contains zigzags Zj and Zj +l . Then the zigzag 
swapping operation (applied to a nowhere zero flow on T) transforms the pair Zj,Zj +l 
into Zj, Zj + %, so the resulting graph T' is not a tree, as it has two edges entering the node 

D. The reduction applied in the proof of part "if" of Theorem ^ in case m > can 
be described in terms of flows. Moreover, the language of flows is convenient to develop a 
more general sort of reduction and to demonstrate some additional properties. To explain 
the idea, consider X G SC(X \ A,0 n , v) and g as in (j2Hj). assuming that A is nonnegative. 
Let P be a path in H beginning at a node (0, s) of the layer L , ending at a node (n,t) 
of the layer L n and such that the minimum a of values of g on the edges of P is nonzero. 
Choose p G Z and a' G 1R. satisfying < s +p < m and < a' < a and change g by moving 
the path P with weight a' at distance \p\, to the right of left depending on the sign of p. 
Formally: define P' to be the path containing the node (i,j + p) for each node (i, j) of P 
and transform g into g' := g — a'x P + a'x P> . This transformation does not change the sum 

1 Note also that for integer points the zigzag swapping operation produces Bender-Knuth's involution, 
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in (|28|). and therefore, the resulting array X' := r )~ l {g') satisfies v x> = v. When p > 
(p < 0), the row derivative <9X' is obtained from dX by increasing (resp. decreasing) by a' 
the entries corresponding to the horizontal edges of the grid G lying between the paths P 
and P'; the tuples X x and A are changed accordingly 

Using such operations, one can transform g more globally, still preserving v: decompose 
g into the sum of path flows ot q x Pq (ot q > 0), q = 1, . . . , N, and move each path P q to the 
left so that the resulting P' q begin at the node (0,0). This gives an array X' with dx'^ = 
for % — 0, . . . , n and j = % + 1, . . . , % + m, i.e., in essense, X' is equivalent to a A-array. One 
can deduce that the first n entries of the tuple A' := X x ' are expressed as follows: 

n+m 

K = Yl lt A *' A *+i] n [Ai,At-fc+i]| for k = l,...,n, (29) 

t=k 

denoting by \[a, b]\ the length b — a of a segment [a, b] and letting Xj := for j > m. Con- 
versely, given A, A, u, define the n-tuple A' by (f2U|) and consider a A-array X' 6 SC(X', n , v). 
Then one can determine a special path decomposition for j(X') and move each path at a 
due distance to the right so as to obtain a flow determining a o- arra y in SC(X \ A, n , i/) 
(moreover, A' is integer when A, A are such and one can maintain flow and array intergality 
under the transformation). This gives a constructive way to reduce the trapezoidal case 
to the triangular one. The tuple A' is weakly decreasing and it just represents the vertex 
generating vector for the permutohedron mentioned in the Introduction. 

Next we explain the idea of deriving Theorem ^ from results in jHJ E] (mentioned in 
Remark 2 in Section We use the equivalence between £A,-arrays of size (n,m) and 
functions on the node set of the corresponding grid G = (V,E). Given tuples A, A, /i, u, 
let us choose a positive integer c and replace //, v by //, v' defined by /i- := //j — ic and 
v[ := Vi — ic, i = 1, ... ,n. This turns the polytope SC(X \ A,/i, v) into SC[X \ A,//, u') 
(each array X in the former polytope corresponds to X' defined by x[j := Xy — 1 ^ ^ c) ; for 
brevity, we denote the latter polytope by C. When c is large enough, C consists of fully 
concave arrays, and we can apply results on the corresponding discrete concave functions. 
The second part of Theorem ^ follows from a result in jH] (in fact, shown there for any 
convex grid) which in our case reads: if A, A, //, v' are integer and if C ^ 0, then C contains 
an integer point. 

The first part of Theorem ^ follows from a combinatorial characterization for the exis- 
tence of a discrete concave function under prescribed boundary data (we use its extension 
to an arbitrary convex grid given in [H]). It uses a notion of puzzle (originally introduced 
for A-grids in [H]). This is a subdivision II of the grid into a set of little triangles and little 
rhombi (the union of two little triangles sharing an edge), along with a 0,1-labeling of the 
edges of G occurring in the boundaries of these pieces, satisfying the following properties: 

(i) for each little triangle r in II, the edges of r are all labeled either by or by 1; 

(ii) for each little rhombus p in II, a side edge of p is labeled 1 if clockwise of an obtuse 
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angle, and if clockwise of an acute angle. 

Then a necessary and sufficient condition on the non-emptiness of C (in £^-case) is that 
each puzzle II satisfies the inequality 

— A( J) + fi'(K) — v'{L) > 0, (30) 

where /, J, K, L are the sets of edges labeled 1 in the lower, upper, left and right sides 
of G, respectively. To show the necessity is rather easy, as follows. Let C ^ and let 
x G C (considering x as a function on V). The discrete concavity of x implies that for 
each little rhombus p with obtuse vertices u,u' and acute vertices v,v', one has q{x,p) : = 
x{u) + x{u') — x(v) — x(v') > 0. When summing up these inequalities for all rhombi in II 
and the equalities (x(t> ) — x{u)) + (x(w) — x(v)) + (x(u) — x(w)) = for all little triangles 
labeled 1, with vertices u, v , w in the anticlockwise order, the terms x(-) for interior vertices 
cancel out and we just obtain (f3T)J) with J, J, K, L to be the sets of edges labeled 1 on the 
corresponding sides. 

When c tends to +oo, the value q(x,p) does so as well (uniformly for all x G C) for 
each little rhombus p, if any, whose smaller diagonal is parallel to the bottom side of G. 
The grow of q{x,p) must cause a similar behavior for the left hand side in (jBTJj) . This 
implies that the puzzles containing at least one of such rhombi p can be excluded from the 
consideration, as they become redundant in verification of the non-emptiness of C. Now 
relation (j3J) in Theorem ^ can be deduced from (jHOj) when the remaining puzzles II are 
considered. 

In conclusion, it should be noted that, using the above reduction to the fully concave 
case and an argument in [2j (where an alternative proof of the integrality theorem from jH] 
is given) , one can show the following sharper version of the last claim in Theorem ^ 

Proposition 2 For integer A, A, p,u, the down hull T> ofSC(\\\,p,v) (i.e., the polyhedron 
SC(X \ A, p, v) — MX ) is integral. 

One can give a direct, relatively simple, proof of this proposition. A sketch: Consider 
a vertex X of T>\ then there is no array X' ^ X in T> with X' > X. Let V\, . . . , be the 
minimal nonempty sets of index pairs such that for q = 1, . . . ,N and for any ij and i'j' 
with i' — i + 1, f G {j, j + 1} and dxij = dxi>j', the set V q contains either both or none of 
ij and i'j'. Let c q := dx^ for ij G V q . Each V q is associated with the corresponding subset 
of horizontal edges in the grid G; let R q denote the union of little triangles containing an 
edge in this subset. Then the interior of each region R q is connected, and each maximal 
horizontal line Li in G (corresponding to the i-th row in dX) intersects R q by a connected, 
possibly empty, set. We say that R q is an intermediate region if it has no edge in the lower 
or upper boundary of G; let for defmiteness R%, . . . , Re be the intermediate regions. One 
shows that if the set V q of nodes of G occurring in the interior of an intermediate region R q is 
nonempty, then one can increase the function a; by a (small) positive constant within the set 
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V q so as to preserve the strip-concavity; the boundary tuples A , A , /x , v are preserved 
automatically. (This relies on the observation that if, e.g., dxtj = dxi-ij and the vertex 
— 1) is in V q , then {i — 1, j — 1) is in V q as well, in view of dxij = dxij-i = dxi-ij-\.) 
Therefore, V q = for all q — 1, . . . in other words, each horizontal line Ci contains at 
most one edge within R q . 

Now associate with R q (1 < q < £) a real variable z q . Let A = (a iq ) be the (n — 1) x £ 
matrix in which ai q is the number of edges of the line L{ occurring in R q . Form the linear 
system Az = b, where for i = 1, . . . , n — 1, hi is equal to Xi^ +m — x i0 minus the sum of values 
dx^ over all ij concerning the edges of non-intermediate regions. Then for the numbers 
c q as above, the tuple z := (ci, . . . , q) is a solution to this system. Note that each bi is 
an integer. (Indeed, each of the above values dx^ is equal to some entry of A or A, which 
is an integer; Xio and Xi^ +m are integers as well.) Also A is a 0,1-matrix and the ones in 
each column go in succession, i.e. A is an interval matrix. So A is totally unimodular 
(cf. [TIH Section 19.4]) and must have full column rank (otherwize Az = has a nonzero 
solution and we can represent X as the half-sum of two other points in SC(X \ A,/i, z/)). 
Then ci, . . . , q are integers, as required. 
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